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1 Introduction 



Among various predictions of QCD, a distinguished one is the existence of the QCD 
phase transition at the temperature T c ~ 150MeV. The phase at the high- 
temperature side, T > T c , is called the quark-gluon plasma (QGP) phase. According 
to the standard big bang scenario, the QGP has lived in the early days of the Universe. 
Efforts of reviving the QGP in the present day are making and it is expected to be 
realized soon. 

Theoretical approaches to the QGP physics fall into four main categories; lattice 
simulation, perturbative approach, effective-theory approach and phenomenological 
approach. These approaches are complementary to each other and each of them has 
its own advantage. 

The purpose of this paper is to give a comprehensive review of the perturba- 
tive approach. [I], |2|, |3|, |], |6|, 0] Thanks to the asymptotic freedom of QCD the 
coupling constant g decreases with temperature T and/or density and then at high 
temperature/density the perturbative approach comes to be an powerful device. 

Thermodynamic properties of a QGP system in thermal and chemical equilibrium 
is characterized by 

Z = Tre-^ 3 -^ (J3 = l/T), 
\nZ = PV/T, 

where Z, H, P and V are, in respective order, the grand-partition function, QCD 
Hamiltonian, pressure and volume of the system, and \i is the chemical potential 
being conjugate to the quark number Q. [For a concise review of the properties of a 
QGP, I refer to Ref. 8).] The thermal average of a quantity Q is defined by 

(n) =Trne-^ H ->" Q) /z. 

Rates of reactions taking place in a QGP are computed through thermal Green func- 
tions, which are defined as the thermal average of the relevant products of field 
operators. 
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2 Perturbative thermal field theory: |t], @, g], |5], 
|6j| Feynman rules 

Traditional approach starts with taking in-fields in vacuum (T = 0) theory as a basis 
of a Fock space. 

2.1 Imaginary- time or Matsubara formalism 

This formalism is convenient for calculating no-leg amplitudes (free energy, grand- 
partition function etc.) and two-point functions. 

• Propagator: 

Po + P 2 + m2 

Here P = (po,p) with p = 2nnT [7r(2n + 1)T — for gluon or FP ghost [quark] 
(n = ■ ■ •, —2, —1, 0, 1, 2, • • •). The form of the "numerator factor" A/"(P) is the same 
as in Euclidean vacuum theory. 

• Vertex: To an iV-particle vertex, is assigned 

V...(27r) 3 i5 n , 5(p), 

where n = Y^f=\ n ii P = Hf=i Pi an d V... stands for the factor that can be read off 
from the interaction Lagrangian. 

• Internal momentum: 




2.2 Real-time formalism 

This formalism allows us to directly compute iV-point functions. The theory is for- 
mulated by introducing a contour C = C\ + C2 + C3 + C4 in a complex time plane: 
U —*tf (Ci), tf —*tf—ia (C 3 ), tf — ia^ti — ia (C 2 ) and ti — ia — ► ti—ift (C 4 ), where 
< u < /3. Then, the limit tj/j — > =1=00 is taken. As far as thermal amplitudes are 
concerned, the contributions from the contour segments C3 and C4 may, in a sense, 
be ignored. H p], Thus the formalism turns out to be a two-component theory; 
the field whose time argument is in the segment G\ (C 2 ), (f>i (02), is called the type-1 
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(type-2) or physical (thermal-ghost) field. Then, in this formalism, propagators, ver- 
tices and self-energy parts enjoy 2x2 matrix structure. Theories with different a's 
constitute an equivalent class of theories. Physical-field amplitudes are independent 
of a. 

• Propagator: 

iM{P)A{P). (2.1) 

Here P = (po, p) with p$ real and [upper (lower) suffix refers to gluon and FP ghost 
(quark)] 

A(P) = M±(P) A F (P) M ± (P) , (2.2) 
A F (P) = diag[A F (P), -A* F (P)] (A F (P) = l/(P 2 -m 2 + *()+)) ,(2.4) 



2V...(2tt) 4 5 4 (P) 



n +^°) = > MPo) = ^R^) + l • (2 - 5) 

iAn, iAi2, «A2i and 1A22 are Fourier transforms with respect to (Re (xq — yo), x — y) 
of, in respective order, (T(p(x)<p(y)) , r (<p(y)(f)(x)), ((p(x)<p(y)) and (T<f)(x)(f)(y)), where 
is the adjoint of 0, r = +(— ) for gluon and FP-ghost (quark) and T is the anti- 
time-ordering symbol. 

• Vertex 

' 1 
v -1 

where P = J2f=i Pi- Note that the vertex matrix is diagonal. The (1,1) component 
is called the type-1 vertex, which is same as the vacuum-theory counterpart, and the 
(2, 2) component is called the type-2 vertex. 

• Internal momentum: 

d 4 P 

(27T)*' 



Notes: 

1) Due to thermal radiative corrections, full gluon- and quark-propagators split 
into several pieces. However their structure remains || essentially to be the same as 
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the corresponding bare propagator, Eqs. ( |2.2| ) - (|2T5| ) (cf. below). Incidentally, (each 
piece of) the self-energy part takes the form 

E(P) = M ± \P) ( S ^ P) J Mg\P) . (2.6) 

2) Thermo field dynamics || f|, (A two-component theory formulated with 
canonical quantization): For each field 0(x) in the original Hamiltonian, introduce 
its "copy" (j)(x). is essentially the field obtained from through time- reversal 
operation, [0 (0) corresponds to 0i (^2) i n the time-path ordered formalism outlined 
above. Note that 02 is in the contour segment C2, in which Re t flies backward, 
+00 — > —00.] Then, a quasiparticle field (fi (i = 1,2) is introduced, in momentum 
space, through 




0(P) = M- 1 (P)0(P) 0(P) = 

<^(P) = $(P) M-\P) , = (0(P), 0(P)) , (2.7) 

where r = +[— ] for gluon and FP ghost [quark]. Note that, in general, Tp i is not the ad- 
joint of ifi. The so-called thermal vacuum is introduced by ^f l \x) \ 0) = (0|^ ^\x) = 
(£1 = +, £2 = — ), where "(+)/(—)" indicates the positive/negative frequency part. 
With this preliminaries, we see that (O|T00|O) = iAfA, Eq. ( |2.1| ). Thus, as far as 
the perturbation scheme is concerned, both thermo field dynamics and the time-path 
formalism summarized above are equivalent. 



3 Reaction-rate formula [|tO| , [tjj 

To avoid inessential complications, I take a heat bath composed of massless scalar 
fields. Reactions taking place in the heat bath are of the following generic type; 

$(Pi) + • • ■ + $(P n ) + heat bath -> $(Q0 + ■ • ■ + $(Q m ) + anything . (3.1) 

Here $ is a nonthermalized heavy scalar particle. [Generalization to other cases is 
st raight forward . ] 
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The reaction rate R reads 
d 

0=1 



V [ n 2 *°rf qj /( 27 r)3 



R 



m J ' 

(3.2) 

where p i0 — Ei — \Jpf + m 2 etc. In Eq. ( |3.2|) , A is an amplitude evaluated in the 
Keldish variant (er = in § 2.2) of real-time formalism for the "process," 

+ . . . + $ l(Pn ) + $ 2 ( Ql ) + ... + $ 2 (g m ) 

-> $ 2 (Pi) + • • ■ + $ 2 (p n ) + $i(g x ) + • • • + $ i(g m ) , 

where the suffix £ i' (i = 1, 2) refers to type-z field. 



Addenda: 

1) A in Eq. ( |3.2[ ) is not an absolute square of some amplitude, in contrast to the 
case of vacuum theory. 

2) When $(-P;) [$((5j)] in the reaction ( |3.1| ) is a thermalized particle, the factor 
n B (Ei) [l + n B (Ej)} is to be multiplied to the right-hand side (rhs) of Eq. (j3^2j). Here 
tib is the Bose distribution function. 

3) The formula ( ^.2| ) is valid [^] even for a finite cube system as far as the periodic 
boundary condition is employed for the single-particle wave function basis. 

4) In the limit T — >• 0, the formula ( |3.2| ) reduces [|ll| to the formula obtained 
through Cutkosky or cutting rules. In particular, for n = 2 and m = 0, the formula 
goes to the optical theorem and, for n = 2 and m — 1, the formula goes to the Mueller 
formula for the corresponding inclusive reaction. 

5) Applying the formula ( |3.2|) to the reaction ( |3.1|) , where <E>'s are constituent 



particles of the heat bath (cf. the second item above), one can derive [0] the detailed- 
balance formula. Namely, the rate ([T^) is equal to the rate for the inverse process to 

m. 



6) Cutting rules: The reaction-rate formula (|3.2j ) is derived from the "first-prin- 
ciple formula" 

R oc Tr e-W-^S'S/Tr e~^ R -^ , (3.3) 
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where S is the S'-matrix element in vacuum theory for the process 

$(Px) + • • • + $(P n ) + {0's} - $(Q0 + ■ • • + $(Q m ) + {0's} . 

Here 0's are constituent particles of the heat bath. 

In what follows, we depart from the scalar theory and keep in mind some general 
theory. The (1,1) component of the thermal propagator A n (P) in A, Eq. ( |3.2| ), 
has the following three roots; (i) the T = propagator connecting the vertices V\ 
and t>2 in S, Eq. (j3lj|), which carry the momentum P from v\ to t> 2 , (ii) a particle 
of momentum p is absorbed from the heat bath into the vertex i> 2 an d a particle of 
the same momentum p is emitted from the vertex v\ into the heat bath and (iii) an 
antiparticle of momentum — p is absorbed from the heat bath into the vertex v\ and 
an antiparticle of the same momentum — p is emitted from the vertex v 2 into the heat 
bath. Both end-point vertices of A n (P) are of type-1, which come from the vertices 
in S in Eq. Q. 

The roots of the (2,2) component of the thermal propagator A 22 (P) in A are 
obtained from above by S — > S*. Both end-point vertices of A 22 (P) are of type-2, 
which come from the vertices in 5* in Eq. ( |3.3|) . 

The (2, 1) component of the thermal propagator A 21 (P) in A has the two roots; 
(i) a particle of momentum p is emitted into the heat bath from a vertex v± in S, Eq. 
( p.3| ), and a particle of the same momentum p is absorbed from the heat bath into 
a vertex t> 2 in S*, Eq. ([H]), and (ii) an antiparticle of momentum — p is absorbed 
from the heat bath into v\ and an antiparticle of the same momentum — p is emitted 
from t> 2 into the heat bath. The vertex v\ (t> 2 ) in S (S*) goes to the type-1 (type-2) 
end-point vertex of A 2 i(P). 

The roots of the (1, 2) component of the thermal propagator A 12 (P) are obtained 
from those for A 21 (P) by S <-> S*. 



This inspection leads us to introduce the thermal cutting rules: |10| Cut all the 
lines A 12 's, A 2 i's, A^P's, and A 2 p's. [The superscript "(T)" refers to the Independ- 
ent part.] 



Through an application of the above cutting rules to some reaction-rate formula, A 
in Eq. ( |3.2j ) is divided into several subparts, each of which corresponds either to S or 
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to S*, Eq. ( |3.3| ), and the interpretation of them in physical terms is straightforward. 

Finally it is worth mentioning that the calculational rules of evaluating absorptive 
part of a generic thermal amplitudes are settled in Ref. 12). Finite-temperature 
generalizations of cutting rules are discussed in Refs. 12) and 13). 



4 Hard-thermal-loop resummation scheme in hot 
QCD jS m 

When formally higher order correction to an (one-particle irreducible) amplitude is of 
the same order of magnitude as the lowest-order counterpart, a resummation of the 
"correction" is necessary. This is the case for classes of amplitudes whose all external 
momenta are soft, P M = O(gT). The relevant diagrams are the one- loop diagrams 
with hard loop momentum, Q\ oop = 0(T), so is named the hard-thermal loop (HTL). 

The computation of 2-point amplitudes or the self-energy parts has been carried 
out long ago. Let us summarize the result. 

Gluon: In a covariant gauge, the full gluon propagator may be decomposed as || 



A' F T (P) = 1/(P 2 -U T (P)) 



(P) - -- 

K } XP 2 + iO+ 

A' F L (P) = 1/(P 2 



-c(P) 



,(4.i; 



P 2 + i0+ 
Hl(P)), (4.2) 



where A is the gauge parameter and Vt [Pl] is the projection operator onto the 
transverse [longitudinal] or chromomagnetic [chromoelectric] sector. The third term 
on the rhs of Eq. ( |4.1|) is the gauge term. Explicit form of Vj!j L , D^ u and C^ u is 
given, e.g., in Ref. 6). 

The HTL contribution reads 



n T (P) 
MP) 

c(P) = 0, 

2 i 

m„ = - 



3 2 
-~m n 
2 9 



PoP 2 
2p 3 

Po 



In 



In 



P 



T 



1 + 



Nf 
6 



Po+P 
\Po-P, 

' po+p \ 
^Po-p) 



7T 2 T 2 



p2 



(4.3) 



where P is soft and Nf is the number of quark flavors. We have assumed the common 
chemical potential // for all Nf (anti)quarks. Note that Ut{P) and Hl(P) are even 
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functions of pq. Observe that II(P)'s are of 0(g 2 T 2 ), the same order of magnitude as 
the bare counterpart P 2 . 



Characteristic features: 



Gl) Landau damping. Im ln[(po +p)/(po — p)} 7^ for space-like P M , P 2 < 0. 

G2) Static limit. Debye screening mass appears in the chromoelectric sector, 
^-l{Po — 0, p) — 3m 2 . On the other hand, no screening mass appears in the chromo- 
magnetic sector, Ht{po = 0, p) = 0. The last fact indicates that, for some amplitudes 
that diverge (in naive perturbative calculation) due to the infrared singularity in the 
chromomagnetic sector, the screening is not sufficient for the amplitudes to converge. 

G3) Dispersion curve [the (positive) solutions, p = uj T / L {p), to P 2 — H T / L (P) = 0]. 
The mode with P 2 — Ul(P) = 0, being absent in vacuum theory, is called the plasmon. 
Po > Trig' The solutions ut/l(p) exist for real p, showing the propagating modes. 
u t/l(p) — P > 0, oj T / L {p = 0) = m g and, for p » m g , uj t {j>) ~ p + 3m 2 /4p and 
^Lip) ~ P + 2p e~ 2p2 / 3m 3 . The group velocities vt/l(p) = d^T/L{p)/dp are positive. 
Po < m g : The solutions po = ut/l(p) exist for pure imaginary p, showing the damping 
modes. 

Quark: The HTL-resummed soft-quark propagator takes the form, 



*S F (P) 



7°-7-p/P + 7° + 7-p/P 



D+(P) 



D-(P) 



m„ 



D±(P) = -(PoTp) + ^ 

2p 

— T 2 + 

6 I 7T 2 T 2 



Po+P 



1 T 5 In 

P \Po-P, 



±2 



(4.4) 
(4.5) 
(4.6) 



Note that D + (—p ,p) = —D_(p ,p). Observe that the HTL contribution, the second 
term on the rhs of Eq. ( |4.5|) , is of O(gT), the same order of magnitude as the bare 
counterpart —po ± p. The 2x2 matrix propagator is related to *Sf{P) through 
M-{P)*S F {P)M-{P) (cf. § 2.2), where *S F {P) = diag[*S F (P), - (*S F (P))*}. Taking 
the complex conjugate, in obtaining the (2,2) component of *S F (P), does not apply 
to the Dirac matrices. 
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Characteristic features: 



Ql) Landau damping as in the case of gluon. 

Q2) Static limit. Debye-like screening mass, D±(p = 0, p) = ±[p + rn 2 q /p\. 

Q3) Dispersion curve [the (positive) solutions, po = uj±(p), to D±(P) = 0]. The 
mode with D_(P) = 0, which is absent in vacuum theory, is called the plasmino. 
Both modes are the propagating modes. u±(p) — p > 0, u±(p = 0) = m q and, for 
p >> m g , uj + (p) ~ p + m 2 q /p and uo-(p) ~ p + 2pe~ 2p2 ^ m ^ . The group velocity of the 
+ mode, v+(p) = du + (p)/dp, is positive. The group velocity V-(p) of the plasmino 
shows an odd behavior. At p — 0, t>_ is negative and, as p increases, V-(p) increases 
across v _ = and approaches u_(p) = 1. Furthermore, at large p, the residue Z_ of 
the plasmino pole damps exponentially, Z_ ~ 2(p 2 /mg)e _2p2//m 9. 

Let me summarize the prominent features of HTL amplitudes. 

1) The HTL contributions, i.e., the contributions that are of the same order of 
magnitude as the lowest-order counterparts, arise in an iV-gluon amplitude (N > 2) 
and a quark-ant iquark-iV-gluon amplitude (N > 0). For an amplitude including 
external FP-ghost lines, which appears in a covariant gauge, the HTL contribution 
does not appear. 

2) The HTL amplitudes are gauge independent. 

3) Kinetic-theory approach leads to the same result. 

4) In vacuum massless QCD, no HTL has arisen. This is because the gauge (chiral) 
invariance of the theory protects a gluon (quark) from getting corrections to the mass. 

5) The free part of the Lagrangian is modified so that the modified one (cf. §5) 
yields the HTL-resummed amplitudes. This means, among others, that, for soft 
modes, the in-field basis in vacuum theory, which is taken as the basis of perturbation 
theory, is not the good basis. It should be noted that the in-fields are irreducible 
representations of Poincare group, which is a symmetry group of vacuum theory. 
However thermal field theory does not enjoy the Poincare symmetry, so that the 
above result is not unnatural at all. In this relation I refer to Ref. 16). 

6) Let G be the exact amplitude with soft external momenta and H be the HTL 
contribution to G. In contrast to the case of vacuum theory, {G — H) / H = O(g). 
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7) HTL iV-point amplitudes satisfy the Ward-Takahashi relation in a stronger 
sense than in vacuum theory. 



From these properties, especially from the last one, one can construct the iV-point 
HTL amplitude (N > 3) from H {2 \ 

5 Effective action [[17], |t§, [tg] 

Having obtained HTL iV-point functions, one can construct an effective action *S, 
which is a generating functional of HTL iV-point amplitudes. *S is the leading con- 
tribution to iSefi defined by 

e i5 cfr = [ V^V$VAe iS , 

./hard modes 

where S is the QCD action. Various forms for *S are available, from which I reproduce 
here the one obtained in Ref. 19): 

q f Y a Y@ r V* 

*<S = —m] J d*xFr{x)( {Y - Ffox)) -m 2 q J d'x^x) (jy^l^)) , 

Y» = (1, k/fc) , 

where D g / q is the covariant derivative acting on the gluon/quark field. 



*S is also deduced [20] from the kinetic-theory approach. 

Various properties of *<S have been disclosed. For interested readers, I refer to the 
literature. 



1) Equation of motion and its solution. [[21 

2) Conserved quantities via Neother's theorem or other means. [p22 

3) Similarity to Chern-Simons theory. |Tj| 23, 24 

4) Classical nature of *S. E4 
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6 Hard modes [26j with \P 2 \ < 0{g 2 T 2 ) 

For illustration of the point, let me cut out the portion from a HTL gluon iV-point 
amplitude with a quark loop, 

S l3 {P + K){-y- l rS ]k {P). (6.1) 

Here Sij is the component of the bare quark propagator constituting the HTL, 
so that P is hard ~ T. K is the momentum of an external gluon and is soft ~ gT. 
Consider, e.g., the sector i — j — 1. Equation (|6.1|) contains 5(P 2 )/(K + P) 2 . When 
the external momentum K is on the mass shell k = ±k, this term develops well- 
known mass singularity at p ■ k = ±pk, 1 / (K + P) 2 oc 1/(1 + p • \t/pk), which, upon 
integration over p, reflects on the logarithmic divergence of the HTL amplitude. This 
is the well-known mass singularity, which appears when the momentum (K + P) M 
can kinematically reach the light cone, (P + K) 2 = 0, on which the bare propagator 
1/(K + P) 2 diverges. This observation leads us to analyze hard propagators near the 
light cone. 

Let us analyze the one-loop self-energy parts, P — > Q + (P — Q) — > P, with P the 
hard external momentum. Recall that, in the case of self-energy part with P soft, 
the hard Q region (HTL) had yielded the dominant contribution. The soft modes 
and the hard modes are "different" modes. By contrast, for hard P, Q and/or P — Q 
are hard. When Q [(P — Q)] is hard, one should use the self-energy-part-resummed 
propagator for the Q (P — Q) line, the self-energy part which we are to evaluate. [As 
a matter of course, when Q [(P — Q)\ is soft, the HTL-resummed effective propagator 
should be used for the Q [(P — Q)) line.] Thus, we are lead to compute the self-energy 
part in a self-consistent manner. 

Here I display the result of the calculation, which is valid to leading order at 
logarithmic accuracy lnl/g >> 1. 

Gluon: The self-energy-part resummed hard-gluon propagators °A^(P)'s in the 
covariant gauge read, with obvious notations, (cf. Eqs. ( f4.1| ) and (f4.2|)): 

♦A£(P) - € M -AjTT^- ( \ ( 6 - 2 ) 

2p p - e(po)\p + 3m 2 g /Ap\ + ze(p )7T 

♦A£(P) ^ (P) ~ j^-^ (6.3) 
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°A£(P) 
A£(P) 

It = 



1 



" A P 2 + iO+ 
-0, 

H9- 1 : 



O 



lnlng 
v ln^" 1 



0{g 2 T) . 



(6.4) 
(6.5) 
(6.6) 



Here A 9 F (P) is the FP-ghost propagator. Above forms are valid in the following 
regions: Im A£(P); ||p | - p\ < O^Tln^" 1 ), Re A T F (P), A L F (P), A 9 F (P), A C F (P); 
0(g 3 T) < ||po| — (p + 3m 2 /4p)| < 0(g 2 T In g" 1 ). 2x2 matrix propagators A's are 
related to A F (P)'s through A = M + A F M + . 

Let us see how the bare propagators are changed through resummation of the self- 
energy part. I take Im °A^(P). The bare form Im A^(P) = — (-n/2p) 5(\po\ —p) turns 
out to be the "smeared" function Im °A^(P), which is peaked at \po\ = p + 3m 2 /4p 
with width j T . Note that m 2 g /p = 0(g 2 T) while = 0(g 2 T In g -1 ), so that » 
m 2 Jp at logarithmic accuracy. For Re °A^(P), similar observation may be made. 



Equations ( |6.2|) - (|6.6|) show that, as in the case of soft modes, the in-field basis 
in vacuum theory is not adequate for the transverse-gluon mode of hard P M with 
P 2 ~ 0. 

Quark (/x = 0): The self-energy-part resummed 2x2 quark propagators reads 

%i{P) * E K*S$\P) (Pr = (l,Tp)) = 1,2) 

Re°S[ T 1 \P) = -Re <> S ( j r 2 \P) 

1 Po ~ 4Po)(P + m 2 q /p) 



Im^S^iP) 



' 2 [p - e(po)(p + m^/p)] 2 + 72 



-7re(po) 



2 ~ ^(^) 



<g(r) 



12/21 ( p ) - -™t(Po) [0(TPo) - 
1 7, 



n F (p)]%(P) 



%(P) 



ir\p - e(po)(p + m 2 /p) + 7 



7 , = £C F T ln^ 1 ) 



2 1 

-i s 



1 + 



lnlng 
, lnc/" 1 



+ 0(/T) 



where the Keldish variant of the real-time contour has been used. 
Similar observation to that in the case of gluon may be made. 
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Substitution of the self-energy-part-resummed propagators screen the above-men- 
tioned mass singularities and renders divergent integral finite. 

7 Application to the computation of physical 
quantities 

Various physical quantities like thermal reaction rates are classified as follows. 

a) Computation in naive perturbation theory yields a finite result. 

b) Naive perturbation theory leads to a diverging result due to infrared (IR) 
singularity which turns out to be finite within the HTL resummation scheme. 

cl) Same as above b) but it still diverges due to the IR singularity. 

c2) Same as above cl) but the divergence is due to the mass singularity. 

Examples of b) are the rate of hard photon and of hard photon-pair productions, 
energy loss of a particle, damping rate of a particle at rest etc. A typical example 
of cl) is the damping rate of moving particle and an example of c2) is the soft- 
photon production rare. A particle at rest feels only chromoelectric field and the 
IR singularity present in the computation within the naive perturbation theory is 
screened by the Debye mass. On the contrary, a moving particle feels chromomagnetic 
field also and, due to the absence of chromomagnetic screening mass, the screening 
at the IR region is not sufficient to render the diverging integral finite (cf. item G2) 
in § 4). 

Hot QCD (rriquark << T) has only two parameters g and T. Then, there arises 
natural hierarchy of scales: T (hard), gT (soft), g 2 T (super soft), ... . 

Naive perturbation theory is valid at the hard region. HTL-resummation scheme 
deals with the soft region. Noting that the quantities classified into a) above receive 
a little contribution from the IR region, one can say that such quantities detect the 
"physics" in the hard region. The quantities classified into b) detect the "physics" at 
the soft as well as hard regions. The quantities belonging to cl) detect the super-soft, 
soft and hard regions. The scheme that deals with super-soft region in a consistent 
manner is not settled yet. In view of the fact that (G — H)/H — 0(g) (cf. item 6) at 
the end of § 4), settlement of this scheme is an urgent issue but there is still a long 
way to go toward the solution. 
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On the other hand, mass-singularity issue seems to be relatively easy to resolve. 

Let me mention, in turn, the damping rate of a moving particle and the soft-photon 
production rate (cf. cl) and c2) above). 

Damping rate: £|/ Within the HTL-resummation scheme, the rate diverges at the 
IR end. It is expected that, at the next-to- leading order, the self-energy part acquires 
screening mass of 0(g 2 T) or "something" which screens the IR singularity. If this is 
the case, the diverging factor ln(gT/0 + ) turns out to be \n(gT/0(g 2 T)) ~ In t7 1 . 

In hot QED, however, it is generally believed that, in any order of perturbation 
series, no magnetic mass is induced. |27| In the IR region, the B loch- Nor dsieck 
approximation, 7^ — > u M (with u M the four velocity), is known to work. Employing 
this approximation scheme, it has been shown |23] that a moving hard electron damps 
according to oc e - aTt M m et) (f. time, m e = eT/3: the QED counterpart of Eq. ( j4.3|) ). 

Soft-photon production rate: The dominant contribution to g^U^ (II^ the pho- 
ton polarization tensor) comes from an one-loop diagram with soft loop momentum. 
|2T| Since all the relevant momenta are soft, one should use HTL-resummed effective 
quark propagators and HTL-resummed photon-quark vertices. As seen at the be- 
ginning of § 6, the HTL photon-quark vertex diverges logarithmically ~ ln(gT/0 + ), 
because the external photon momentum is on the mass shell. 

According to the general argument in §6, substitution of self-energy-part-resum- 
med hard-quark propagators *S"s for the bare ones make diverging result finite, 
ln(<7T/0+) - HgT/0(g 2 T)) ~ ln^ 1 . 

However, the above substitution violates the Ward-Takahashi relation, which in- 
dicates that there must be important vertex corrections. Lebedev and Smilga have 



shown [|30| that the corresponding diagrams are the (resummation of) ladder dia- 
grams. This yields the additional contribution to the soft-photon production rate, 
which coincides with the above contribution to leading order at logarithmic accuracy. 



31] Whether or not this analysis can be generalized to a generic reaction rate or 



thermal amplitude that belongs to the category c2) above is an open question. 
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8 Beyond the hard-thermal-loop resummation 
scheme 



As mentioned in the last section (cf. also the item 6) at the end of § 4), toward 
establishing a next-to-leading-order resummation scheme is still a long way. Here I 
simply enumerate, without comment, some of the work made toward this end. 



1) Next-to-leading order computation of (chromoelectric) Debye mass. |32| 

2) Next-to-leading order computation of plasmon frequency. |33| 

3) Next-to-leading order computation of the gluon vacuum polarization tensor 
IV (P) with soft P. 

4) Next-to-leading order correction to the dispersion laws (cf. items G3) and Q3) 
in §4). |H 

5) Self-consistent determination of chromomagnetic mass. |3lJ 

6) Improved effective action. p7 | 

Although not directly related to the subject of this review, I enumerate the fol- 
lowing important achievements in the field of thermal field theory. 



1) Hot QED and hot scalar QED. pTI, |38 



2) Calculational scheme of grand partition function or pressure. |39| 
Finally I mention the extensions to the nonequilibrium thermal field theory. For 
dealing with systems that are quasiuniform near equilibrium or quasistationary, tra- 
ditional approach uses |5], the Keldish variant of real-time formalism. As far as the 
computation of reaction rates are concerned, almost all the machineries of equilibrium 
thermal field theory hold as they are. An important one that does not hold is Eq. 
(ETBD, which causes the appearance of pinch singularity in self-energy-part-inserted 
propagators. [41 Two approaches are devoted to this issue. 1) It has been shown 



| 42]| that such singular contributions can be resummed (see also Ref. 5)). Application 
of this result to the hard-photon production rate is made. |^3j 2) A renormalization 
theory constructed through renormalizing number densities is proposed. MM This 
theory is same in structure as the equilibrium thermal field theory, so that no pinch 
singularity appears. 

Thermo field dynamics as mentioned in §2 (cf. item 2) after Eq. (pl]|)) is gener- 
alized to the nonequilibrium case. 0] Recall that the Bogoliubov matrix M± in Eq. 
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Q2.7| ) defines the quasiparticle fields <£>'s. Now M± and then the quasiparticle picture 
depend on space-time coordinates. From Eq. ( |2.7| ), we see that the time derivative 
of y?'s receives two contributions, one coming from the time derivative of the original 
fields 0's is governed by the Hamiltonian and the other comes from the time derivative 
of M±, through which the thermal energy is introduced. Then, through renormal- 
ization procedure of propagators, time dependence of M± or the number density 
is determined. The determining equation turns out to be a generalized Boltzmann 
equation. 

9 Conclusion 

The structure of perturbative hot QCD is far more complicated than the perturbative 
vacuum QCD. Naive perturbation scheme, which is formulated using in-field basis 
(in vacuum theory) in a Fock space, is valid only at the low level of the "QCD 
mountain" (the short wave-length or the hard region A ~ 1/T). At the high level of 
the mountain (the long wave-length or the soft region A ~ 1/gT), naive perturbation 
scheme breaks down, which means that the in-fields are not the good basis at this 
level. The perturbation scheme that works here is the HTL-resummation scheme. 
Again this scheme does not apply at yet higher level of the mountain (the longer 
wave-length region). Continuous efforts aiming at establishing the new resummation 
scheme that works at this level are making. The goal is, however, still far a way. 

There are "ravines" along the light cone (mass singularities) here and there in the 
QCD mountain. Techniques of pass over these ravines are not completely settled yet. 

Comprehensive analysis of rates of various reaction taking place in nonequilibrium 
system, as well as the development of the theoretical framework of nonequilibrium 
quantum-field theory per se, have begun. 
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